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toluene-benzene  dimers  using  a  combination  of  atom-atom  and  multipolar  potent¬ 
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a  molecular^quadrupole-quadrupole  interaction  term  gives  the  reported  parallel 
displaced  C^^. benzene  dimer  geometry.  The  excited  state  behavior  of  the 
benzene  dimer  can  be  qualitatively  understood  based  on  this  calculation.  The 
toluene-benzene  dimer  is  well  modeled  bv  these  terms  in  addition  to  a  dipole-^ 
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20.  qiiadrupole  interaction  term.  Both  parallel  displaced  and  perpendicular 

configurations  are  obtained  in  agreement  with  the  experimental  data.  Using 
exponential-six,  dipole-dipole,  quadrupole-quadrupole ,  and  dipole-quadrupole 
interactions  the  toluene  dimer  can  be  satisfactor ilv  modeled.  The  calcula¬ 
tions  generally  suggest  that  excitation  of  parallel  displaced  dimers  can  lead 
to  rearrangement  (and  possible  excimer  formation)  in  the  excited  state. 
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I.  INTRODUCTION 

One  of  the  most  important  properties  of  a  van  der  Vaals  (vdW) 

cluster  is  its  geometry:^  cluster  geometry  expresses  not  only  the 

symmetry  and  steric  properties  of  a  system,  but  also  its  energtics. 

Such  observables  as  spectral  shifts  and  binding  energy  are  sensitive 

2 

to  changes  in  the  geometry  of  a  cluster.  Since  the  equilibrium 
geometry  of  a  cluster  is  the  result  of  a  delicate  balance  of  a  variety 
of  intermolecular  interactions,  a  knowledge  of  cluster  geometry  can 
lead  to  a  better  understanding  of  the  relative  sizes  of  the  various 
contributions  to  the  overall  cluster  interaction  energy. 

Computer  modeling  of  the  intermolecular  energy  of  a  vdW  cluster 
has  been  shown  to  be  a  very  useful  contribution  to  the  process  of 

3 

deducing  cluster  geometry  from  spectroscopic  data.  Minimum  energy 


configurations  of  the  clusters  of  interest  are  obtained  for  a  number 
of  different  intermolecular  potentials.  A  geometry  consistent  with 
all  spectroscopic  observables  is  accepted  as  the  sufficient  one.  In 
addition,  this  process  suggests  a  form  for  the  intermolecular 
interactions. 

Spectroscopic  data  for  aromatic  dimers  ((benzene)^,  (toluene)2, 

tol-ben  are  not  always  sufficient  to  determine  the  geometry  of  the 

dimers  uniquely.^  For  example,  a  parallel-displaced  (pd)  dimer  is 

.4  &! 

suggested  by  the  2-color  MS  data  for  (ben;^  »  less  information  ; 

■d 

is  obtained  from  the  spectra  of  (tol)^  and  tol-ben.  The 


spectroscopic  data,  in  addition,  do  not  shed  light  on  such  questions 


as,  why  do  dimers  have  certain  geometries,  why  are  some  dimer  red 


_ _  •iJ'lity  Codes 
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Bhifts  large  and  others  small,  why  do  dimer  excited  state  lifetimes 
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differ  from  monomer  lifetimes,  or  how  can  excimers  form  upon 
excitation  of  ground  state  dimers.  The  computer  modeling  of  aromatic 
dimers  has  enabled  us  to  address  these  and  a  number  of  other  questions 
discussed  in  this  pap-r. 

Recent  reports  from  this  laboratory  have  demonstrated  the 
usefulness  of  potential  energy  calculations  for  benzene-solvent, 
toluene-solvent,  and  n-propyl  benzene-solvent  clusters.'  An  atom- 
atom  exponential-six  (exp-6)^  form  for  the  intermolecular  potential 
provides  reasonable  geometries  for  the  ben(S)^,  tol(S)^,  and 
npb(S)^  (x<3)  systems.  However,  the  aromatic  dimers  are  not  well 
represented  by  an  exp-6  form  of  the  intermolecular  potential:  the 
predicted  geometries  are  in  general  not  consistent  with  the 
spectroscopic  data.  The  approach  taken  by  Williams  to  improve  the 
intermolecular  pair  potential  involves  placing  partial  charges  on  the 
atoms. ^  Although  reasonably  successful,  this  approach  is  not  ideal 
because  some  clearly  incorrect  low  energy  structures  are  obtained  and 
the  structures  are  quite  sensitive  to  the  size  of  the  partial  charges 
chosen.  We  have  instead  employed  the  exp-6  form  supplemented  by 
molecule-molecule  multipolar  terms  such  as  quadrupole-quadrupole 
(V^^),  quadrupole-dipole  dipole-dipole 

interactions  as  appropriate.  While  this  procedure  is  subject  to 
limitations  as  well,  it  at  least  has  the  advantage  that  experimental 
values  of  many  of  these  moments  are  available  in  the  literature.  By 
exploring  the  relative  importance  of  such  terms  and  by  estimating  the 


magnitudes  of  the  multipolar  interactions  in  the  excited  state,  dimer 
geometries  and  energetics  can  be  elucidated. 

This  paper  reports  the  results  of  computer  modeling  of  {ben)^, 
(tol)^,  and  tol-ben.  Spectroscopic  data  for  these  aromatic  dimers 
demand  the  use  of  forms  of  the  intermolecular  potential  more 
complicated  than  exp-6  atom-atom  potentials.  The  choice  of  which 
extra  terms  to  include  is  governed  by  symmetry  and  the  observed  data. 
II.  PROCEDURES 

The  computer  modeling  algorithm  used  to  generate  the  aromatic 

dimer  energy  and  geometry  is  similar  to  methods  previously  employed 
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for  solute-solvent  clusters.  However,  in  this  study  five  different 
expressions  for  the  intermolecular  energy  are  used  (see  Table  l).  The 
first  involves  only  an  exp-6  potential  with  parameters  determined  by 
Williams:^  these  terms  are  present  in  all  other  forms  of  the 
potential.  The  parameters  are  chosen  to  reproduce  the  heat  of 
simulation  of  the  benzene  crystal.  Since  the  interactions  in  the 
crystal  and  the  aromatic  dimers  are  in  general  similar,  such  parameters 
should  serve  as  a  reasonable  first  approximation  of  the  dimer  inter¬ 
action.  The  second  energy  expression  is  the  improved  potential  of 
Williams^  which  places  partial  charges  on  the  atoms.  The  remaining 
three  expressions  incorporate  molecule-molecule  terms  appropriate  for 
(ben)^,  tol-ben,  and  (tol)^  as  a  replacement  of  the  coulomb  part  of 
the  intermolecular  potential.  Expressions  for  ^(jq* 

7 

V, ,  are  derived  from  the  text  of  Hirschfelder,  Curtis  and  Byrd 
dd 

0 

utilizing  the  quadrupole  moment  of  benzene  and  dipole  moment  of 


toluene.  The  numerical  values  of  the  coefficient  for  the  exp-6, 

^qq’  ^dq  ^dd  sane  in  all  expressions.  The 

quadrupole  moment  of  toluene  is  assumed  to  be  the  same  as  that  of 
benzene.  The  geometry  of  benzene  is  that  used  by  Williams^,  and  the 
geometry  of  toluene  is  taken  from  the  crystal  stmicture.^^  The 
center  of  the  multipole  moments  of  toluene  is  displaced  from  the 
center  of  the  aromatic  ring  (.464A)  towards  the  methyl  group  at 
approximately  the  center  of  charge. 

In  order  to  estimate  the  dimer  spectral  shifts  and  obtain  some 

idea  about  excited  state  energetics  and  structure,  excited  state  dimers 

are  also  modeled.  The  general  approach  employed  involves  changing  the 

multipole  moments  of  the  molecule  excited:  the  exp-6  parameter  could 

not  be  readily  varied  in  a  systematic  fashion.  The  excited  state  of 

toluene  is  modeled  in  two  separate  calculations:  one  calculation 

doubles  both  the  dipole  and  quadrupole  moments  of  toluene  and  the  other 

increases  the  dipole  moment  by  a  factor  of  8.  This  latter  rather  large 

1 2 

change  is  suggested  by  a  recent  calculation.  The  actual  change  in 
the  moments  of  toluene  most  likely  falls  between  these  two  limits. 

The  success  of  these  intermolecular  potential  calculations  is 
largely  due  to  the  separability  of  the  interactions.  It  has  been 

1 3 

shown  that  spectral  shift  scales  with  polarizability  and  dipole 
moment. Since  to  a  first  approximation  the  exp-6  terms  model  the 
dispersion  forces  only,  the  multipolar  terms  can  be  added  as 
appropriate  without  re-optimizing  the  exp-6  parameters. 


^V-V 


The  routine  for  minimizing  the  intermoleculer  energy  of  the 

3 

dimers  is  the  sane  as  that  described  previously  for  clusters.  The 


energy  for  a  particular  configuration  is  calculated,  then  the  ligand 
is  moved  a  step  in  one  of  the  six  (3  translational  and  3  rotational) 
degrees  of  freedom  and  the  energy  is  recalculated.  If  the  resulting 
energy  is  lower,  the  step  is  repeated;  if  the  energy  is  not  lower  the 
ligand  is  stepped  back  and  another  degree  of  freedom  is  explored.  A 
minimum  energy  configuration  can  thus  be  achieved.  Care  must  be  taken 
to  start  the  cluster  in  a  number  of  sufficiently  different  , 

configurations  to  ensure  that  all  of  the  local  minima  in  the  potential 
surface  (configuration  space)  are  located.  Atomic  positions  for  all 
of  the  calculated  geometries  are  available  upon  request.  Only 
pertinent  data  for  the  geometries  presented  in  figures  are  given  in 
the  text. 

III.  RESULTS  ARC  DISCUSSION 

A.  Benzene 

Spectroscopic  Data-  The  spectroscopic  data  for  (ben)„ 


reveals  much  about  the  geometry  of  the  dimer.  A  0^  transition  is 
clearly  evident  for  the  dimer,  whereas  this  transition  is  symmetry 


forbidden  in  the  monomer.  The  dimer  must,  therefore,  have  less  than 
C^  symmetry.  Since  only  one  peak  is  observed  in  the  origin  region, 
only  one  type  of  benzene  molecule  in  the  dimer  is  present  in  appreciable 
amounts;  both  benzene  molecules  in  the  dimer  must  be  spectroscopically 


equivalent  (see  the  discussion  below  for  (tol)^  and  tol-ben).  A  C, 
parallel-displaced  (pd)  geometry  is  consistent  with  the  data.^ 


V  > 
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The  observed  spectral  shift  for  (ben)^  is  -42  cm”  •  Comparison 
of  this  spectral  shift  with  previous  results  allows  a  qualitative 
estimate  of  the  relative  displacement  of  the  benzene  molecules  in  the 
dimer.  For  clusters  studied  previously  (ben(S)^,  tol(S)^,  npb(S)^)  , 
the  spectral  shift  is  found  to  be  proportional  to  the  polarizability  of 
the  solvent  and  the  displacement  of  the  solvent  away  from  the  center  of 
the  aromatic  ring.  The  large  polarizability  of  the  benzene  "solvent" 
implies  that  the  relative  displacement  should  be  large  ('u2A)  to  produce 
such  a  relatively  small  (-42  cm  spectral  shift. 

Computer  Modeling-  The  benzene  dimer  is  modeled  using  three 
different  forms  for  the  intermolecular  potential  (labeled  I,  II,  III 
in  Table  I).  The  exp-6  potential  calculation  predicts  only  one 
geometry  for  (ben)^  (see  fig.  l):  the  two  molecules  are  centerea 
over  one  another  in  a  staggered  configuration.  This  geometry  of  the 
dimer  preserves  the  6-fold  symmetry  axis  of  benzene  and  consequently 
must  be  rejected. 

The  second  expression  used  to  calculate  the  geometry  of  (ben)^ 

contains  a  partial  charge  atom-atom  term  (q.q./R.  .).  The  q. 's 

t  j  1  j  1 

have  been  determined  to  fit  the  free  energy  of  the  benzene  crystal^ 

(see  Table  I).  This  expression  predicts  two  stable  geometries  for 
(ben)2  (see  fig.  1);  a  "herring-bone"  and  a  pd  geometry.  The  two 
benzene  molecules  in  the  herring-bone  configuration  are  physically 
distinguishable  and  this  geometry  is  therefore  rejected.  The  pd 
geometry  is  consistent  with  all  spectroscopic  data  and  has  a  binding 


energy  within  15  cm”  of  that  of  the  herring-bone  configuration. 


V  • .  ." .'  I  ^'1 V  ■.■'  •?> 


Although  this  expression  for  the  intencolecular  energy  predicts  a 
configuration  consistent  with  the  (ben)^  spectroscopic  data,  it  has 
at  least  two  shortcomings.  First,  the  lowest  energy  configuration 
generated  by  this  calculation  is  not  consistent  with  the  (ben)^ 
spectrum.  The  herring-bone  geometry  should  also  be  observed 
spectroscopically  if  this  second  calculation  is  valid.  Second,  the 
geometries  predicted  by  calculations  which  include  the 
term  are  very  sensitive  to  the  values  chosen  for  q.  and  q.. 

Reducing  the  values  of  q.  and  q,  to  those  which  yield  the  correct 
molecular  quadrupole  moment  eliminates  the  herring-bone  configuration 
and  greatly  reduces  the  ring  displacement  in  the  pd  configuration.  A 
displacement  this  small  would  be  expected  to  produce  a  spectral  shift 
much  larger  than  -A2  cm~^  .  On  the  other  hand,  increasing  the 
can  eventually  generate  a  perpendicular  configuration. 

The  third  potential  (Table  I)  employed  by  modeling  this  dimer 

consists  of  the  exp-6  terms  in  addition  to  a  molecular  quadrupole- 

quadrupole  interaction  This  interaction  results  in  a 

geometry  for  the  (ben)^  completely  consistent  with  the  spectroscopic 

data.  The  term  in  the  potential  incorporates  the  correct  value 

for  the  benzene  quadrupole  moment  and  causes  the  two  benzene  molecules 

o 

to  be  displaced  in  the  X  (C-H  bond)  direction  by  2.5A  (fig.l).  A 
displacement  this  large  is  consistent  with  a  -42  cm”^  spectral  shift. 
The  calculated  geometry  is  pd  and  has  symmetry.  The  two  benzene 
molecules  are  indistinguishable,  only  one  stable  configuration  is 
predicted,  and  the  spectra  are  well  modeled  in  all  respects. 
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Excimer  Formation-  Potential  calculations  for  the  benzene 
dimer  shed  some  light  on  the  formation  of  the  suggested  (ben)^ 
excimer.  Spectroscopic  data  point  clearly  to  the  formation  of  an 
excimer:  the  lifetime  of  all  features  in  the  spectrum  falls  from  100 

ns  in  the  isolated  molecule  to  40  ns  while  the  fluorescence  quantum 
yield  falls  significantly  (probably  a  factor  of  10  ).  The  expected 
geometry  for  the  excimer  is  or  parallel  undisplaced  and 

closely  spaced  (<3-5A  ring  separation).^  The  above  calculations 
present  a  possible  mechanism  for  the  formation  of  the  excimer  from  the 
pd  ground  state  dimer.  The  calculations  just  described  (potentials  I 
and  III,  Table  I)  suggest  that  the  displacement  of  the  two  benzene 
molecules  in  the  benzene  dimer  is  due  essentially  to  the  V  term. 

qq 

Excitation  of  benzene  to  the  state  causes  the  molecular 

quadrupole  moment  to  vanish  and  thus  the  V  term  contribution  to 

qq 

the  binding  energy  is  "zero"  for  the  excited  state  benzene  dimer.  The 
dimer  can  then  collapse  to  the  parallel  undisplaced  configuration 
usually  associated  with  excimer  formation  without  a  significant 
barrier. 

The  suggested  dimer-excimer  transition  still  presents  somewhat  of 
a  puzzle:  the  excited  state  lifetime  of  the  dimer  is  not  consistent 
with  the  very  low  quantum  yield  for  emission  and  the  diner  absorption 
features  are  sharp  (  2.5  cm”^ )  with  large  Franck-Condon  overlap  for 
the  low  lying  vdV  modes.  Such  sprectral  features  are  not  characteristic 
of  absorption  to  high  lying  levels  of  an  excimer  well.  The  present 
calculations  may  suggest  a  possible  mechanism  for  this  behavior.  The 
pd  ground  state  can  undergo  a  transition  to  a  pd  excited  state 


stabilized  by  small  induced  quadrupole  moments  in  benzene.  If  this 
pd  excited  state  lives  for  'v>2ps  and  then  collapses  to  a  parallel 

5 

undisplaced  configuration,  the  experimental  data  can  be  rationalized. 
B.  Toluene-Benzene 

4 

Spectroscopic  Data-  The  spectroscopic  data  for  tol-ben 
reveal  two  sets  of  peaks  in  the  region  of  the  toluene  0^  transition 
(tol  -ben  Oq).  One  set  of  peaks  has  a  small  spectral  red  shift 

(<25  cm~^),  a  low  intensity,  and  is  composed  of  at  least  two 

*  0 

separate  origins.  The  other  set  of  peaks  in  the  tol  -ben  0^  region 
has  a  large  spectral  shift  (-250  cm~^),  only  one  obvious  origin 
followed  by  a  progression  of  vdW  stretching  and  bending  modes,  and  is 

about  10  times  more  intense  than  the  peaks  with  the  small  shift. 

*  1 

The  spectrum  in  the  tol-ben  6^  region  shows  only  one  set  of 
features.  Presumably  the  lower  overall  intensity  of  the  spectrum  in 
this  region  precludes  observation  of  the  less  shifted  set  of  peaks. 

The  geometry  described  as  exhibiting  a  large  spectral  shift  for  the 

*  0  *  • 

tol  -ben  0^  region  has  a  much  smaller  shift  in  the  tol-ben  region 

( -40  cm  ^ ) . 

Computer  Modeling-  The  computer  modeling  of  tol-ben 
using  the  exp-6  form  for  the  potential  energy  predicts  three  stable 
geometries:  two  perpendicular  ( J_ )  configurations  and  one  pd 
configuration  (see  fig  2).  The  spectroscopic  characteristics  expected 
for  the  _[^  configurations  are  those  of  the  less  shifted  set  of  features 

*  A  , 

in  the  tol  -  ben  0^  region.  The  spectral  shifts  for  the  calculated 
geometries  should  be  small,  assuming  the  general  relationship  between 
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spectral  shift  and  overlap  of  the  "solvent"  with  the  aromatic  syster 
of  the  solute  is  valid.  The  prediction  of  more  than  one  J_  geometry 
is  consistent  with  the  observation  of  more  than  one  origin  in  the  less 
shifted  set  of  features.  Finally,  the  weak  intensity  expected  for  the 
less  shifted  peaks  correlates  well  with  the  observed  intensities:  the  J_ 
geometries  are  less  stable  than  the  pd  geometry  (see  figure  2). 

The  tol-ben  pd  geometry  calculated  with  the  exp-6  potential  (fig. 

2)  is  not  entirely  consistent  with  the  observed  spectrum.  Although  the 
large  binding  energy  of  the  pd  configuration  serves  as  a  rationalization 
for  the  high  intensity  of  the  features  at  -250  cm”^  from  the  toluene 
Oq,  the  observed  shifts  of  tol-ben  transitions  are  not  consistent  with 
the  pd  geometry  predicted  by  the  exp-6  calculation.  One  would  anticipate 
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red  shifts  for  the  tol-ben  transitions  of  much  more  than  40  cm”  and 

*  -1 

red  shifts  for  tol  -ben  much  less  than  250  cm  :  in  fact,  expecta¬ 
tions  based  on  this  potential  form  are  for  equal  shifts  for  the  excitation 
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of  either  molecule  in  this  geometry. 


Addition  of  and  terms  to  the  exp-6  potential  form 
(Table  I,  potential  IV)  provides  calculated  geometries  consistent  with 
the  observed  binding  energies,  spectral  shifts,  and  intensities.  The  J_ 


geometries  generated  by  this  improved  potential  are  similar  to  those 

calc’  ''.ated  with  the  exp-6  terms  only.  However,  the  pd  configuration 

generated  with  potentials  I  and  IV  are  quite  different  (fig.  2).  The 

greater  displacement  of  the  two  molecules  in  the  dimer  obtained  with 

the  exp-6  +  V  +  V,  reduces  the  r-system  overlap  between  benzene 
^  qq  dq 

*  1 

and  toluene  such  that  the  expected  tol-ben  6^  red  shift  relative 
to  benzene  6^  is  more  like  that  found  for  the  (ben)^  6q. 


-1  * 
Furthenaore,  the  -250  cm  spectral  shift  for  the  tol  -ber. 

Oq  can  now  be  readily  understood.  The  dipole  moment  (end  presumably 

the  quadrupole  moment)  of  toluene  is  expected  to  increase  upon 

1 2 

excitation;  perhaps  by  as  much  as  a  factor  of  eignt.  This 
increase  in  multipole  moment  increases  the  binding  energy  of  the 
complex,  and  therefore  the  8p<*ctral  shift.  The  trends  can  be 
approximately  modeled  by  increasing  the  mult  ole  interaction  terms 
until  the  desired  change  in  intemol ecul ar  potential  energy  is 
achieved.  Increasing  the  quadrupole  and  dipole  moments  of  toluene  by 
a  factor  of  two,  or  increasing  the  dipole  moment  by  a  factor  of  eight 
produces  an  increase  in  the  binding  energy  of  200  cm”'.  The 
equilibrium  geometry  of  this  model  "excited  state"  complex  is  very 
similar  to  the  ground  state  geometry  and  therefore  is  consistent  with 
the  large  Franck-Condon  overlap  for  the  first  few  vibronic  levels. 

Thus  a  modest  increase  in  the  quadimpole  moment  of  toluene  (benzene) 

and  a  reasonable  increase  in  the  dipole  moment  upon  excitation  would 

,  _1  .  * 
account  for  the  large  (-250  cm  )  spectral  shift  for  the  tol  -ben 

0^  transition.  The  actual  change  in  multipole  moments  upon 

excitation  is  probably  within  the  limits  expressed  above  (i.e.,  Q"  < 

Q’  <  2Q"  and  D"  <  D'  <  8D"). 

Calculations  of  the  tol-ben  geometry  using  potential  IV  (Table  l) 

yield  an  additional  stable,  nearly-pd  (10°  from  parallel)  configuration 

o 

in  which  the  benzene  is  displaced  2.6A  towards  the  methyl  group  of 
toluene.  Mo  spectroscopic  features  in  the  tol-ben  spectrum  can  at 


present  be  assigned  to  this  dimer  geometry.  Since  the  term  is 


positive  in  this  configuration,  an  excited  state  rearrangenent  to  the 


nore  stable  pd  geometry  discussed  above  could  well  occur.  Absorption 

* 

to  this  state  (tol  -ben)  would  then  be  weak,  broad  and  blue  shifted 

with  respect  to  the  tol  0^;  it  could  be  buried  in  the  region  between 

♦  0  *  0 

the  tol  -ben  0^  and  the  tol-ben  0^.  It  is  also  possible  that 
this  configuration  may  undergo  thermal  rearrangement  in  the  ground  state 
to  the  more  stable  geometry  depicted  in  figure  2,  as  the  calculated 
barrier  to  rearrangement  is  estimated  to  be  '^50  cm  \ 

Excicer  Formation-  The  lifetime  data  for  the  peaks 
in  the  tol-ben  spectrum  show  that  an  excimer  forms  upon  excitation  of 
some  of  the  tol-ben  vibronic  levels.  The  data  are  consistent  with  a 
barrier  to  excimer  formation  of  approximately  900  cm~^ ;  roughly  900 
cm  of  toluene  excess  vibrational  energy  must  be  excited  before 

4 

appreciable  excimer  formation  occurs. 

The  formation  of  the  tol-ben  excimer  upon  excitation  of  benzene 

can  be  explained  as  before  (see  Benzene  Dicer  -  Excimer  Formation). 

The  quadxnipole  moment  of  benzene  vanishes  in  the  excited 

state  and  the  V  and  V,  terms  drop  out  of  the  intermolecular 
qq  dq 

potential  energy  expression  leaving  only  the  exp-6  form  (Table  l). 


Without  and  terms,  the  minimum  energy  geometry  is  only 
slightly  displaced  and  excimer  formation  should  readily  occur  (figure 


On  the  other  hand,  both  the  spectroscopic  data  and  the 

calculations  (potential  expression  IV,  Table  l)  indicate  that  the  pd 

*  0 

geometry  of  tol-ben  changes  very  little  at  the  tol  -ben  0^  level. 

• 

Therefore,  the  pd  geometry  of  tol  -ben  should  feature  the  center  of 

o 

the  benzene  ring  shifted  about  2A  from  the  center  of  the  toluene 


aromatic  ring.  This  displaced  configuration  is  stabilized  by  the 
multipolar  terms  in  the  energy  expression  and  a  barrier  to  excimer 
formation  is  thereby  quite  reasonable.  Conceivably,  nearly  1000 
cm”^  of  excited  state  molecular  vibrational  energy  could  be  required 
to  transfer  enough  energy  to  the  vd¥  modes  to  bring  the  benzene  into  a 
position  favorable  for  excimer  formation. 

C.  Toluene  Dimer 

Spectroscopic  data-  Little  is  revealed  about  the  geometry 

0  4 

of  the  toluene  dimer  from  its  spectrum  in  the  0^  region.  Two  sets 
of  absorption  features  can  be  identified.  One  set  exhibits  a  small 
spectral  shift  (-10  cm~^ )  and  has  relatively  sharp  features.  The 
progression  in  the  vdW  stretch  for  this  set  of  features  shows  a  good 
Franck-Condon  overlap  near  the  dimer  0^.  The  other  set  of  peaks 
consists  of  a  broad,  structured  background  extending  far  to  the  red. 
Pew  assignments  can  be  made  in  this  region  of  the  spectrum  and  several 
geometries  are  likely  to  contribute  to  the  congestion. 

Computer  Modeling-  Computer  calculations  with  potential  V 
(Table  I)  yield  a  number  of  different  geometries  for  the  toluene 
dimer,  pd,  J.,  and  even  nearly  -  J.g®ometries  appear  to  be  stable 
because  of  the  interplay  between  the  various  orienting  terms  in  the 
full  potential.  In  addition,  the  nearly  -  J_ configurations  are  highly 
dependent  on  the  values  of  the  multipolar  moment  parameters  and  may 
not  actually  be  stable  configurations.  Below  we  will  discuss  only 
the  more  certain  pd  and  configurations  of  (tol)^.  Figure  5  contains 
the  calculated  geometries  for  the  pd  and  J_  configurations. 


The  configuration  displayed  in  figure  3  is  stable  for 


reasonable  values  of  the  molecular  moments.  The  two  toluene  molecules 
are  distinguishable  in  the  J_  configuration  and  should  give  different 
spectral  shifts:  the  horizontal  (fig.  3)  toluene  would  be  expected  to 
yield  the  smaller  spectral  shifts.  Consequently,  the  sharp  features 
in  the  {tol)^  Oq  spectrum  are  assigned  to  excitation  of  the 
horizontal  (fig  3-)  toluene  in  the  X  configuration.  Excitation  of  the 
vertical  toluene  should  lead  to  features  with  larger  spectral  shifts, 
presumably  obscured  by  the  broad,  more  red-shifted  features  mentioned 
above. 

The  ground  state  pd  geometry  predicted  by  an  exp-6  potential  is 

very  different  from  that  predicted  by  the  full  potential  (exp-6  ♦ 

V  +  V ,  +  V  ,  +  V,.),  as  shown  in  figure  3.  The  geometry 

qq  dq  qa  dd  ^ 

calculated  with  this  full  potential  is  very  sensitive  to  the  assumed 
value  of  the  quadrupole  moment  of  toluene.  Clearly,  the  spectral 
shift  associated  with  the  pd  geometry  is  larger  than  that  of  theX 
geometry.  Doubling  of  the  toluene  dipole  moment  causes  substantial 
reorientation  of  the  toluene  dimer:  doubling  the  quadrupole  moment 
has  little  effect  on  this  geometry.  Increasing  the  dipole  moment  of 
toluene  by  a  factor  of  8  produces  a  qualitatively  similar  geometry. 
This  calculated  dramatic  change  in  equilibrium  configuration  upon 
excitation,  and  concomitant  increase  in  binding  energy,  can  serve  to 
rationalize  the  diffuse  nature  of  the  more  red  shifted  features  in  the 
(tol)^  spectrum.  The  ground  state  dimer  may  have  a  maximum  Franck- 
Condon  overlap  with  a  dense,  high  energy  portion  of  the  vdW  vibrational 
manifold  of  the  excited  state:  a  broad  relatively  unstructured 
spectrum  can  thus  be  generated. 
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Sxcimer  Formation-  The  lifetime  data  for  some  of  the  toluene 
dimer  0^  features  suggest  excimer  formation.  The  sharp  features 
shifted  by  -10  cm  ^  from  the  toluene  0^  do  not  show  lifetime 
shortening,  but  the  broad  more  red  shifted  features  do.  The  computer 
modeling  of  the  toluene  dimer  geometry  is  consistent  with  this 
observation;  excimer  formation  is  most  likely  to  occur  for  a  pd 
rather  than  a  J. geometry.  The  calculations  further  indicate  that  a 
barrier  to  excimer  formation  might  occur.  If  absorption  is  to  high 
energy  vdW  modes  (as  discussed  previously),  the  excess  vdW  energy 
could  lead  to  excimer  formation  dispite  the  barrier. 

IV.  CONCLUSIONS 


An  atom-atom  exp-6  potential  is  not  adequate  for  the  computer 


modeling  of  the  aromatic  dimers  (ben)^,  (tol)^*  and  tol-ben.  The 


addition  of  a  quadrupole-quadrupole  interaction  term  (V  )  for  the 


modeling  of  (ben)„  results  in  a  pd  geometry  consistent  with  the 


spectroscopic  data.  The  magnitude  of  the  spectral  shift  (-45  cm”^) 
is  consistent  with  the  predicted  displacement  of  the  two  rings. 
Moreover,  since  the  displacement  of  benzene  in  the  dimer  is  due  to  the 


V  term  and  the  quadrupole  moment  vanishes  in  the  excited  state, 

qq 


the  (ben)„  upon  excitation  moves  to  an  undisplaced  parallel  geometry 


which  allows  excimer  formation. 


The  tol-ben  dimer  can  be  modeled  by  the  addition  of  V  and 

qq 


terms  to  the  basic  exp-6  potential  form.  TheJ.  geometries 


computed  within  this  format  should  produce  features  consistent  in 


*  0  0 
shift  and  intensity  with  the  tol  -ben  0^  features  near  the  tol  Oq 


transition.  The  large  spectral  shift  of  the  other  set  of  tol  -ben 
0^  features  can  be  associated  with  a  pd  geometry.  The  large  change 


in  binding  energy  upon  excitation  in  the  pd  configuration  is  explained 

by  a  significant  increase  of  the  molecular  dipole  noment  and  a 

moderate  increase  in  the  quadrupole  moment  of  toluene  in  the  excited 

state.  The  large  relative  displacement  of  the  two  rings  in  this  dimer 

even  in  the  excited  state  nay  give  rise  to  the  observed  barrier  to 

excimer  formation.  Excitation  of  benzene  in  the  pd  geometry  of  tol-ben 

causes  the  V  and  V.  potential  terms  to  vanish, 
qq  dq 

The  (tol)2  spectrum  is  qualitatively  consistent  with  the 
geometries  predicted  by  the  calculations  if  ®nd 

terms  are  added  to  the  intermolecular  potential.  The  two  toluene 
molecules  are  distinguishable  in  the  J_  {tol)^  geometry.  Excitation  of 
the  horizontal  toluene  (fig.  3)  in  thej.  geometry  produces  the  less 
red  shifted  features  in  the  0^  spectinim.  The  broad,  more  red 
shifted  features  can  be  associated  with  other  J.  forms  and  the  pd 
configuration.  In  the  pd  configuration  the  increase  in  the  toluene 
dipole  moment  upon  excitation  may  cause  significant  reorientations  of 
the  molecules.  The  absorption  spectrum  of  the  pd  (tol)^  may  then  be 
broad  because  the  Eranck-Condon  maximum  for  the  transition  occurs  for 
a  dense  manifold  of  higher  energy  excited  state  vdV  modes.  Excimer 
formation  is  expected  for  the  pd  but  not  thej.  geometry. 

The  interplay  of  model  calculations  and  spectroscopy  has  proven 
essential  to  the  understanding  of  aromatic  dimers.  Without  the 
calculations  the  spectroscopic  results  are  quite  puzzling.  The  major 
test  of  the  calculations,  of  course,  is  that  they  are  consistent  with 
the  experimental  observations.  Taken  together,  the  calculations  and 
the  data  reveal  a  good  deal  about  the  Investigated  dimers. 
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FIGURE  CAPTIONS 


FIGURE  1 

Minimum  energy  configurations  calculated  for  (ben)^  using  three 
different  expressions  for  the  intermolecular  potential  (Expressions  I, 
II  and  III  in  Table  I).  The  total  intermolecular  potential  energy  (in 
cm”^ )  for  each  of  the  calculated  dimers  is  given  in  the  figure-  The 
exp-6  +  geometry  adequately  models  the  ground  state  dimer  while 
the  exp-6  configuration  models  the  equilibrium  geometry  of  the  excited 
state  before  excimer  formation  (see  text). 

FIGURE  2 

Important  minimum  energy  geometries  for  tol-ben  calculated  using  two 
different  intermolecular  potentials  (expressions  I  and  IV  in  Table  l). 
The  exp-6  +  V  +  V,  potential  shows  three  ground  state 

qq  aq 

geometries  assignable  to  features  in  the  tol-ben  spectrum.  Excited 
state  geometries  are  not  shown.  Binding  energies  are  given  in  cm~^. 
FIGURE  3 

Representative  minimum  energy  geometries  of  (tol)^.  Several  J.  and 
nearly-  configurations  are  found;  pictured  are  the  geometries  for 
which  spectroscopic  features  are  discemable.  The  pd  configuration 
undergoes  substantial  change  in  geometry  upon  excitation.  The  excited 
state  is  modeled  here  by  doubling  the  ^qqf  ^dq’  %d  ^dd 
terms.  Binding  energies  are  given  inc  cm 
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